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Abstract. We give an alternative proof of a result by N. Gantert, Y. Hu and 
Z. Shi on the asymptotic behavior of the survival probability of the branching 
random walk killed below a linear boundary, in the special case of deterministic 
binary branching and bounded random walk steps. Connections with the Brunet- 
Derrida theory of stochastic fronts are discussed. 



1. Introduction 

Consider a real-valued branching random walk, with deterministic binary branch- 
ing and i.i.d. random steps. One possible way of constructing such a branching 
random walk is to start with a rooted infinite binary tree T and a collection (£(e)) e 
of i.i.d. copies of a random variable where the index e runs over the edges of T. 
Given the value of the branching random walk at every vertex u of T is then 

defined as 

S(u):=x+ Yl C(e), 

ed[root,u\ 

where [roof, it]] denotes the set of edges of T lying between u and the root. With 
this definition, for each infinite ray root =: xq, x%, . . . in T, the sequence S(xi)i = Q j i ; ... 
formed by the successive values of the random walk along the ray, is a real-valued 
random walk started at S(root) = x. In the sequel, we use the notation P x to reflect 
the choice of x as the initial value, so that Po is used for the case where the random 
walk starts at zero. 

For v E M, define A^fa) to be the event that there exists an infinite ray 
root =: xq, x\, . . . in T such that S(xi) > vi for all i > 0. A dynamical inter- 
pretation of this event is as follows. Assume that the tree T depicts the genealogical 
structure of a population whose individuals are identified with the vertices, so that 
the population starts at generation zero with a single ancestor represented by the 
root. For each vertex u £ T, S(u) is interpreted at the location on the real line of 
the individual represented by u. Now assume that, in every generation n > 0, any 
individual whose location is strictly below vn is killed, and removed from the popu- 
lation together with all its progeny. In this setting, the event A^fa) corresponds to 
the fact that the population survives the killing process forever , i.e. that the n— th 
generation is non-empty for every n > 0. 

The authors would like to thank J. Quastel for suggesting that the approach developed in [27, 28 
might be exploited to re-derive the asymptotic behavior of the survival probability of the branching 
random walk. 
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It turns our that, under suitable regularity assumptions on the distribution of 
C, there exists a critical value v* such that ¥q(A OCi (v)) > for all v < v* , while 
IPo(^4oo(^)) = for all v > v* (see e.g. [9]). In this paper, we are interested in how 
fast ¥q(A 00 (v)) goes to zero when v < v* goes to v*. To state a precise result, we 
first introduce some notations and assumptions on the distribution of £. Our first 
and rather drastic assumption is that ( has bounded support, or, in other words, 
that there exists a bounded interval [£_, £ + ] such that 



is well-defined and finite for all i £ R, and is C°° as a function of the parameter 
t. Our next (and last) assumption on the distribution of C is thafl there exists a 
t* s]0, oo[ such that 



and a recent result by N. Gantert, Y. Hu and Z. Shi [21] gives the following descrip- 
tion of the asymptotic behavior of the survival probability as v goes to v*: 

Theorem 1. For v < v* , one has the following asymptotic behavior as v —> v* : 



In fact, the theorem quoted above is less general than the one proved in [2~T] . 
since we have assumed deterministic binary branching and bounded random walk 
steps, whereas [21] allows for a more general branching mechanism, and possibly 
unbounded steps provided that the Laplace transform is finite in a neighborhood 
of zero. In this paper, we shall prove the following slightly improved version of 
Theorem [U 

Theorem 2. For v < v* , one has the following asymptotic behavior as v —> v* : 



The main novelty in the present paper lies in the method of proof, which is 
completely different from that developed in [21]. Indeed, the proof given in |21j is 
based, among other things, on a first-second moment argument, using a change-of- 
measure technique combined with refined "small deviations" estimates for random 
walk paths, exploiting some ideas developed in [24] in the context of branching 
Brownian motion. On the other hand, our proof relies on the characterization of the 
survival probability of the branching random walk as the solution of a non-linear 

^The meaning of this assumption is discussed in more detail in the Appendix. 




(2) A(t*)-t*A'(t*) = -log(2). 
With these assumptions, the critical value v* is given by 

(3) v* = A'(t*). 
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convolution equation. This idea was in fact used by B. Derrida and D. Simon in [17\ 
[30] to derive a heuristic justification of Theorems [Hand El treating the corresponding 
equation by the non-rigorous methods developed by Brunet and Derrida to study 
stochastic front propagation models. The rigorous treatment given here was inspired 
by the work of C. Mueller, L. Mytnik and J. Quastel \27\ I28j . who deal with a 
continuous-time version of the non-linear convolution equation, as an intermediate 
step in their proof of the Brunet-Derrida velocity shift for the stochastic F-KPP 
equation. Both our argument and the one given in |27[l28j use the idea of comparing 
the solutions of the original non-linear equation to solutions of suitably adjusted 
linear approximations of it. An important difference, however, is that the equation 
appearing in |27U28j is a second-order non-linear o.d.e., for which specific techniques 
(such as phase-plane analysis) can be applied, while such tools are not available in 
our discrete-time setting, so that we had to find a different way of implementing the 
comparison idea. 

The rest of the paper is organized as follows. In Section [21 we introduce the 
non-linear convolution equation characterizing the survival probability, and show 
how, given super- and sub-solutions of this equation, one can obtain upper and 
lower bounds on the survival probability. Section [3] is devoted to the study of 
a linearized version of the convolution equation, for which explicit solutions are 
available when v is close to v*. In Section [J] we explain how the explicit solutions 
of the linearized equation derived in the previous section, can be used to build 
super- and sub-solutions to the original non-linear convolution equation. Finally, 
Section [5] puts together the arguments needed to prove Theorem [21 Section [6] then 
discusses the connection of Theorem [2] and its proof with the Brunet-Derrida theory 
of stochastic fronts and related rigorous mathematical results. In the appendix, we 
discuss the meaning of the assumption that there exists a t* such that ([2]) holds, 
and the asymptotic behavior of the survival probability when this assumption is not 
met. 

2. Equations characterizing the survival probability 

In this section, we explain how to the survival probability can be characterized 
by non- linear convolution equations, and how super- and sub-solutions to these 
equations provide control upon this probability. 

2.1. Statement of the results. Throughout this section, v denotes a real number 
such that 

V < V* . 

Let V : [0, 1] -»■ [0, 1] be defined by 

V>(s) = 2s- s 2 . 
Let H denote the following set of maps: 

% := {h : M — > [0, 1], h non-decreasing, h = on ] — oo,0[ }, 
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and, for all h E H, define T[h) E W. b30 

T(h)(x) := ij){E{h{x + ( - v))), x>0 



T(h){x) := 0, x < 

For all n > 0, let A n (u) denote the event that there exists a ray of length n 
root =: xo,xi, . . . ,x n in T such that S(xi) > vi for all i E [0, nj. Finally, for all 
xGK, let 

q n (x) :=F x (A n (v)), goo (x) :=P«(Ax,(«))- 
Note that, for obvious reasona^, g n and goo belong to %, and go = l[o,+oo)- Given 
hi, h,2 E we say that /ii < /12 when < /i2(x) for all igl. 

The following proposition gives the non-linear convolution equation satisfied by 
q n and goo, on which our analysis is based. 

Proposition 1. For all n > 0, g n +i = T(q n ), and T(q oa ) = q^. 

Proof. Analysis of the first step performed by the walk. □ 

A key property is the fact that being a non-trivial fixed point of T uniquely 
characterizes g^o among the elements of H, as stated in the following proposition. 

Proposition 2. Any r E % such that r ^ and T(r) = r is such that r = q^. 

Our strategy for estimating qoo(x) is based on the following two comparison re- 
sults. 

Proposition 3. Assume that h E % is such that h(0) > and T(h) < h. Then 
q<h. 

Proposition 4. Assume that h E H is such that T(h) > h. Then q > h. 

2.2. Proofs. Let us first record the following simple but crucial monotonicity prop- 
erty of T. 

Proposition 5. Ifhi,h2 E % satisfy hi < /i2, then T(hi) < T(h2). 

Proof. Immediate. □ 

As a first useful consequence of Proposition El we can prove Proposition [H 

Proof of Proposition [|} Since T{h) > h, we can iteratively apply Proposition [5] to 
prove that for all n > 0, T n+1 (h) > T n (h), whence the inequality T n {h) > h. On 
the other hand, since h E H, we have that h < l[o,+oo[> whence the inequality 
T n (h) < T n (lr +oo [). Since T n (l[ 0j +oo[) = Qn by Proposition [Q we deduce that 
q n > h. By dominated convergence, lim n _> +00 q n (x) = qoo(x), so we finally deduce 
that h < goo- □ 

We now collect some elementary lemmas that are used in the subsequent proofs. 

Lemma 1. One has that P(( < v*) < 1. 



2 The fact that T(h) £ H is a straightforward consequence of the fact that h £ H and that ifi is 
non-decreasing. 

For instance by an immediate coupling argument. 
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Proof. By assumption, there exists t* > such that A'(t*) = v* . But A'(t*) = 
~w^t w c^ > so that P(Ce 4 *^) = E(v*e t *^). If moreover P(£ < v*) = 1, we deduce from 
the previous identity that P(£ = t>*) = 1, so that A(t*) = which contradicts 
the assumption that A(f) - = - log(2). □ 

Lemma 2. One /ms i/iat goo(0) > 0. 

Proof. In Section[4j for v* — v small enough, we exhibit c_ G % such that T(c_) > c_ 
and c_(x) > for x > 0. We deduce (with Lemma [T]) that T(c_)(0) > 0. Now, 
letting /i := T(c_), Proposition [5] shows that T(h) > h. Proposition 2] then yields 
that goo(0) > h(Q) > 0. This conclusion is valid for small enough v* — v. Since 
clearljQ </oo(0) is non-increasing with respect to v, the conclusion of the lemma is in 
fact valid for all v < v*. □ 

Lemma 3. There exists a constant k > depending only on the distribution of £, 
such that, for all v < v* , goo(0) > k<Zoo(1)- 

Proof. Using the fact that ip(s) > s for all x £ [0, 1], we have that, for all x G [0, +oo[, 
> E(q 00 (x + £ — w)). From Lemma [H we can find r/ > such that P(( > 
v* +rj) > 0. Using the fact that q is non-decreasing, we obtain that qoo(x) > P(( > 
v * + rj)q 00 (x + rj). Iterating, we see that, for all n > 0, qooO) > P(( > v * + r n) n q 00 {nrj). 
Choosing n large enough so that nrj > 1, we get the desired result. □ 

We now prove Proposition [2j 

Proof of Proposition Let r € T~L satisfy T{r) = r and r ^ 0. According to Propo- 
sition HI we already have that r < q^. Our next step is to show that r(0) > 0. To 
this end, let 

D :={x£ [0,+oo[; r(x) > 0}. 
Since we assume that r ^ 0, D is non-empty, and since, moreover, r is non- 
decreasing, D has to be an interval, unbounded to the right. Since v < v* , we 
know from Lemma [T] that P(£ < u) < 1, whence the existence of rj > such that 
P(C — v > rj) > 0. Let x be such that x + r\ belongs to D. Then, x + Q — v belongs 
to D with positive probability, so that E(r(x + £ — t> )) > 0. If, moreover, x > 0, we 
have that: 

r(x) = T{r)(x) = ip(E(r(x + £ — f))), 
so that r(x) > since ^(s) > for all < s < 1. We have therefore proved that 

(-D-ry)n [0,+oo[c L>. 

Since D is a subinterval of [0, +oo[, unbounded to the right, this implies that D = 
[0, +oo[, whence r(0) > 0. Now let 

F := {A > 0;Vx G [0,+oo[, r(x) > Agoo(x)}. 

Since goo(0) > by Lemma [2] and r < q^, F must have the form [0, Ao] for some 
Ao £ [0, 1], and, we need to prove that indeed Ao = 1 to finish our argument. We 
first show that Aq > 0. Since r is non-decreasing, we have that, for all x G M, r > 



4 For instance by coupling. 
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r(0)l[o i + 00 [, whence r > r(0)qoc since l[o,+oo[ > Qoo- As a consequence, Ao > r(0), 
and we have seen that r(0) > 0. Now, using the fact that T(r) = r, Proposition [5l 
and the definition of Ao, we see that 

r = T(r) >T{\ oqoo ). 

For x > 0, we also have that qoo(x) > 9oo(0) > 0, and E{q QO {x + ( — v)) > 0. Since 
Ao > 0, we can write: 

r(-r\ > \ n a ('A T ( X 0Q)( X ) — \ n „ (t\ T ( X 0<1°°)( x ) 
T{X) > MqcoiX) Xoqoo(x) - ^OQcoiX) XoT(qZ)(x) ' 

whence the inequality 

r(x) > X qoo(x) x{E(qoc{x + (-v))), 
where \ is the map defined, for s G]0, 1], by 

^(A s) 2 - \ s 



X(s) ■-- 



AoV'O) 



with the extension x(0) := 1. S ince (/oo is non-decreasing, this is also the case of the 
map x t-t E(q 00 (x + ( — v)). Moreover, % too is non-decreasing on [0, 1], so we get 
that: 

r{x) > \ Q x{E{qoc{C, - v)))q 00 (x). 

Therefore, ^ox(E(,Qoo(( — v))) is an element of the set F. If Ao < 1, the fact that 
E(q 00 (( — v)) > and strict monotonicity of \ show that x{E{q{C,— v))) > x(0) = 1. 
We would thus have the existence of an element in F strictly greater than Ao, 
a contradiction. We thus conclude that Ao equals 1. Therefore one must have 
r > qoo- □ 

Proof of Proposition^ Since T{h) < h, we deduce from Proposition [5] that, for all 
if!, the sequence (T n (h)(x)) n >o is non-increasing. We deduce the existence of a 
map T°°(h) in U such that, for all x G R, T°°(h)(x) = lim n ^ +00 T n (h)(x), and it 
is easily checked by dominated convergence that T(T°°(h))(x) = T°°(h)(x) for all 
while T°°{h){x) < h(x). It remains to check that T°°(h) ^ to obtain the 
result, since Proposition [2] will then prove that T°°{h) = q^. Using the fact that 
h is non-decreasing, we have that, for all x > 0, h(x) > /i(0)lr 0j+oo r(a;). Moreover, 
it is easily checked that, for all A, s G [0,1], tp(Xs) > Xip(s). Using Proposition [5l 
we thus obtain that T n (h)(x) > h(0)T n {l [0:+c<) [)(x), whence T°°(/ i )(x) > h(0) qoo (x) 
by letting n — > +oo. We deduce that T°°(h)(x) j£ since we have assumed that 
h(0) > 0. □ 

Remark 1. The proof of Proposition^ given above relies solely on analytical argu- 
ments. It is in fact possible to prove a slightly different version of Proposition @ - 
which is sufficient to establish Propositions^ and\Q - using a probabilistic argument 
based on the interpretation of the operator T in terms of branching random walks. 
We thought it preferable to give a purely analytical proof here, since our overall proof 
strategy for Theorem is at its core an analytical approach. 
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3. Solving a linearized equation 

According to Proposition^ the survival probability qoc can be characterized as the 
unique non-trivial solution (in the space H) of the non-linear convolution equation, 
valid for x > 0: 

(5) r(x) = ip(E(r(x + C-v))). 

Linearizing the above equation around the trivial solution r = 0, using the fact that 
ip(s) = 2s + o(s) as s — > 0, yields the following linear convolution equation: 

(6) r(x) =2E(r(x + C-v)). 

As explained in Section [5] below, explicit solutions to a slightly generalized version 
of flSJ) are precisely what we use to build super- and sub-solutions to the original 
non-linear convolution equation ([5]). To be specific, the linear convolution equation 
we consider are of the form 

(7) c(x) = e~ a E(c(x + C - x€R, 

where c : R — > C is a measurable map, v is close to v* and e~ a is close to 2. Looking 
for solutions of (J7|) of the form 



(8) c(x) = e 



')..{: 



where <j> G C, we see that a necessary and sufficient condition for (jSj) to yield a 
solution is that 

(9) E (e^-^ = e a . 

Due to our initial assumption on £, we have that 

E (e l *^- v *A = 1/2, 



so that one can hope to find solutions to (jUJ) by performing a perturbative analysis. 
This is precisely what is done in Section 13.11 Then, in Section 13.21 some of the 
properties of the corresponding solutions of ((7J) are studied. 



3.1. Existence of exponential solutions. Consider the extension of the Laplace 
transform of £ when the parameter t € C. Since £ has bounded support, t \— > 
E(exp(t£)) defines a holomorphic map from C to C, and, since E(exp(t*C)) ^ 
we may extend the definition of A to an open neighborhood U of t* in C, by setting 
A(t) := log-E7(exp(i£)) for t & U, using the principal determination of the logarithm. 
We thus obtain a holomorphic map on U. 

Proposition 6. Let a be a holomorphic function defined on a neighborhood of zero, 
such that a(0) = — log(2) and a'(0) = 0. There exists eo > and a map 4> ■ [0, eo] ~~ * 
U such that, for all v G [v* — eq,v*] ; the following identity holds: 

(10) E (e<t>(v*-v)(£-v)^ =e a{v*-v)^ 
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and such that, as e — > 0, 



/ 1t*f 

(11) ^ e) = t * + ^A^F) +0{e) - 

Proof. Given v in the vicinity of v* , we are looking for a t £ U such that 

(12) E (e t{<: - v A = e a{v *- v) . 



For t close enough to t*, and v close enough to v*, we can use the logarithm and 
observe that the equation 

(13) A(t)-tv = a(v* -v). 

is sufficient for (I12p to hold. Expanding A for t £ £7, we have that 

A(t) = A(0 + {t - t*)A'(f) + (t- t*) 2 g(t), 

where g is holomorphic on U and satisfies g(t*) = A"(t*)/2. Plugging ([2]) and ([3]), 
we can rewrite the above expansion as 

A(t) = -log(2) + tv* + (t-t*) 2 g(t). 

On the other hand, using our assumptions on a, we can write 

a(v* - v) = - log(2) + (v* - v)b(v* - v), 

where b is an holomorphic function in a neighborhood of zero such that 6(0) = 0. 
Finally, (fl~3|) reads 

(14) {t . t *f g{t) = {v _ vl{t _ b{v *_ v)) _ 



Observe that, since g(t*) = A"(t*)/2 ^ R_, we can define y/g(t) for t close to 
t* , using the principal determination of the logarithm and the definition y/z = 
exp(log z/2). We can similarly define yt — b(z) for (t,z) close to (t*,0). Now, for 
(t,u) G C X C in the vicinity of (t*,0), consider the equation . 

(15) (t-t*)y/g(f) = uy / t-b(-u 2 ). 

Clearly, if (jl5h holds with u = i^/v* — v when v* and v are real numbers such that 
v < v*, then (|14p holds. Now consider the map H defined in the neighborhood of 
(t*,0) in C x C by 



E(t, u) := (t - f)vf)- u^t-bi-u 2 ). 

Observe that 5(t*,0) = 0, and that the (holomorphic) derivative of H with respect 
to t at (t*,0) is equal to y/g(t*) = v / A"(t*)/2 ^ 0. Identifying C with RxR, 
we can thus view 5 as a smooth map defined on an open set of M 4 , and apply the 
implicit function theorem to deduce the existence of a smooth map / defined on a 
neighborhood of in C such that /(0) = t* and such that, for all u near zero, 

(16) S(/(«),u) = 0. 

We now set <j>(e) = f(iy/e), which yields (fTOj) . Then, one obtains (fTT|) by computing 
the derivative of / at zero from (|16p in the usual way. □ 



SURVIVAL PROBABILITY OF THE BRANCHING RANDOM WALK 



9 



3.2. Properties of the exponential solutions. Now let a, cq and (ft be given as 
in Proposition [H Throughout the sequel, we use the notation 

e = v* — v. 

An immediate consequence of the proposition is that, for all v 6 [v* — €q,v*], the 
map defined on K by x \— > e^ e ' x , solves the equation 

(17) c(x) =e- a(e) E{c{x + (-v)). 

If a(e) 6 M. when e£K (this will be the case in all the examples we consider below), 

then the map x h-> e^ e ' x is also a solution of (fT7|) . where z denotes the conjugate 
complex number of z. Let us set a(e) := 3?(^>(e)) and (3(e) := S(0(e)). Thus, we 
obtain a solution of (fT7|) if we set 

(18) d(x) := e a(e)a: sin(/3(e)x). 

Consider e small enough so that a(e) > and (3(e) > 0. Note that 

d(0) = d(n/(3(e)) = 0, 

and that one has 

d < on [-vr//3(e),0], 
d > on [0,vr//3(e)], 
d < on [Tr/(3(e),2TT/(3(e)]. 

The derivative of d is given by 

d'(x) = a(e)e a ^ x sm(/3(e)x) + (3(e)e a ^ x cos(/3(e)x). 

One thus checks that d attains a unique maximum on the interval [0, ir/(3(e)], at a 
value x = L(e) satisfying 

(19) tan(/3(e)L(e)) = -/3(e)/a(e), 

and that d is increasing on the interval [0, L(e)}. As e goes to 0, we know from (jlip 
that 



/ 1f*e 

(20) a(c) = t* + 0(6), /3(e) = J + 0(e), 

and we then deduce from P^|) that, as e goes to 0, 

(21) L(e) = ^-l/f + o(l), 
whence 



A"(t*) 

(22) L(c) =,^-^ + 0(1). 
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4. Building a super-solution and a sub-solution to the original 

equation 

In this section, we explain how to transform the explicit solutions of the linear 
equation (|7|) obtained in the previous section, into super- and sub- solutions of the 
non-linear convolution equation of Section [2J In the sequel, A denotes a real number 
(whose existence is guaranteed by the assumption that ( has bounded support) such 
that, for all v in the vicinity of v* , 

P(\(-v\ < A) = 1. 

4.1. The super-solution. Let us choose the function a in Proposition [6] as the 
constant function a(e) := — log(2). Then consider the function d(-) defined in (|18p . 
and note that any function of the form 

c(x) := A(e)d(x + A), x€R 

is a solution of ([7]). Now let A(e) be implicitly defined by the requirement that 

c(L(e)-2A) = 1. 

This last condition rewrites more explicitly as 

(23) A(e)e a ^ L ^-^ sin(/3(e)(L(e) - A)) = 1. 

From (|2ip . we have that, as e — > 0, 

sin(/3(e)(L(e)-A))~/3(6)(l/t* + A). 

Combining with ([23]), (J22J) and (J20D, we deduce that, as e -> 0, A(e) > and 



(24) logA(e) = -^^^- + O(loge). 

For notational convenience, we introduce 

C(e) :=L(e) - 2A. 

The next proposition summarizes the properties of c that we shall use in the sequel. 

Proposition 7. For small enough e, the following properties hold. 

(i) c(x) > for allxe [-A,0]/ 

(ii) < c(x) < 1 for all x £ [0, C(e)]; 

(iii) c{x) > 1 /or all x £ [C(e), C(e) + A]. 

(iv) c is non- decreasing on [0, C(e)] 
Moreover, as e goes to zero, 



/ \"(t*)t* 

(25) logc(l) = -Tr^^-^ + 0(loge). 

Proof, (i), (ii), (iii) and (iv) are rather direct consequences of the definition and of 
the analysis of Section [372J As for (|25p . it is readily derived from ()24|) and (|20p . □ 

We now define the map c+ : K — > [0, 1] by 

• c + (x) = for all x < 0; 

• c+(x) = c(x) for all x G [0, L(e) - 2A]; 
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• c+(a?) = 1 for all x > L(e) - 2A. 
Now let V>+ be defined on [0, 1] by 

i/j+(s) := min(2s, 1). 
Proposition 8. For all s G [0, 1], ifi+(s) > ip(s). 

Proof. Immediate. □ 

Proposition 9. For small enough e, one has that c + G % and, for all x > 0, 

c+(x) > if> + (E{c+(x + C-v))). 

Proof. The fact that c + G H is guaranteed by the definition and properties (i) to 
(iv) of Proposition Consider first the case x G [0, C(e)]. By definition, one has 
that 

c(s) = 2E(c(x + (- v)). 
Since, by construction, c + < c on [—A, C(e) + A], we deduce that 

c(x) > 2E{c + (x + ( -«)). 

Since tp+(s) < 2s for all s G [0, 1], we deduce that 

c(x) > ^+(E(c+{x + C-v))), 
whence, remembering that c+{x) = c(x), 

c + (x)>^ + (E(c + (x + (-v))). 

Now, for x > C(e), we have that c + (x) = 1, so that c + (x) > f° r au s ^ [0, 1]- 

In particular, 

c+(x) > ^+{E(c+{x + C - v))). 

□ 

Corollary 1. For small enough e, one has that c+ > T(c+). 

Proof. Direct consequence of Propositions [8] and [9l □ 

4.2. The sub-solution. In this section, we choose the function a in Proposition [6] 
as the function a(e) := — log(2) + e 2 . For notational convenience, we also introduce 

7 (e) :=e-^), h(e) := 7 (e)(2 - 7 (e)). 

As above, we use the function (i(-) defined in p8h . and let 

c(x) := A(e)d(x) = A(e)e a(e):r sin(/3(e)x), 

where -A(e) is implicitly defined by the requirement that 

c(L(e)) = h{e). 

The same argument as in the previous section then shows that, as e — > 0, A(e) > 
and 



(26) logA(e) = -K^j^&L + 0(loge). 

Here are the properties of c that we shall use in the sequel. 
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Proposition 10. For small enough e, the following properties hold. 

(i) c(x) < for all x G [-A, 0]; 

(ii) < c(x) < h(e) for all x G [0,L(e)]; 

(iii) c(x) < h(e) for all x G [L(e),L(e) + A]; 

(iv) c is non- decreasing on [0, L(e)]. 

Moreover, as e goes to zero, 



(27) logc(l) = + O(loge). 

Proof, (i), (ii), (iii) and (iv) are rather direct consequences of the definition and of 
the analysis of Section [372J As for (|27p . it is readily derived from (j26j) and ([20p . □ 

We define the map c_ : R — > [0, 1] by 

• c_(x) = for all x < 0; 

• c_(x) = c(x) for all x G [0,L(e)]; 

• c_(x) = /i(e) for all x > L(e). 
Now define ip_ on [0, 1] by 

ip-(s) := min(j(e)s,h(e)). 

Proposition 11. For small enough e, for all s G [0, 1], ip-(s) < tf)(s). 

Proof. Immediate. □ 

Proposition 12. For small enough e, one has c_ G 7i and, for all x > 0, 

c_(x) < ip-(E(c-(x + C - «)))• 

Proof. The fact that c_ G % is a direct consequence of the definition and of proper- 
ties (i) to (iv) of Proposition [TOl Consider first the case x G [0, L(e)]. By construc- 
tion, one has that 

c(x) = r y{e)E{c{x + Q — v)). 
Since, by construction, c_ > c on [— A,L(e) + A], we deduce that 

(28) c{x) <j{e)E(c-{x + (-v)). 
Since we have c(x) < h(e) and c_(x) = c(x), we deduce that 

c-{x) < ip-(E(c-(x + C - v))). 

Consider now x > L{e). Since c_ is non-decreasing, so is y \- > ^{e)E(c-{y + C — v)). 
We deduce that 

-y(e)E(c-(x + ( - v)) > 7 (e)S(c_(L(e) + C - «)). 
Using again the fact that c_ > c on [— A,L(e) + A], we have that 

7 (e)E(c_(L(e) + (-v))> j(e)E(c(L(e) + ( - v)) = c(L(e)) = h(e). 
Since c_ < h(e), we finally deduce that 

c_(x) < ^_(JS(c_(x + C - «)))• 

□ 
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Corollary 2. For small enough e, one has that c_ < T(c_). 

Proof. Direct consequence of Propositions [IT] and H2J □ 

5. Proof of Theorem [2] 

We now put together the different pieces leading to the proof of Theorem [2j For 
small enough e, Corollary [JJ yields a super-solution of the non-linear convolution 
equation, i.e. a function c + £ rl such that c + > T(c+). By Proposition O we 
automatically have that goo < c+. Now we have goo(0) < 900(1) < c+(l), so that the 
asymptotic behavior stated in (I25p yields that 



loggoo(O) < -J^^ + 0(log(v* - „)). 

Similarly, Corollary [2] yields a sub-solution, i.e. a function c- £ rl such that c_ < 
T(c_), and, by Proposition H] we have that goo > c_. Using Lemma O we obtain 
that goo(0) > k<?oo(1) > K c_(l), so that the asymptotic behavior stated in (]2"7|) yields 
that 

loggoo(O) > -^0^ + OQog(v* - v)). 
This concludes the proof of Theorem [2j 

Remark 2. One may wonder whether it is possible to extend our method of proof 
to a more general situation such as the one treated in [21] , i.e. a fairly general 
supercritical stochastic branching mechanism, and random walk steps whose Laplace 
transform is finite in a neighborhood of zero. We expect our proof to be robust to more 
general branching mechanisms, since the key properties of the resulting non-linear 
convolution equation (monotonicity, qualitative features of the linearized equation) 
should be preserved. On the other hand, it might not be easy to replace the bounded 
support assumption by a less stringent one, since one would have to control how 
large unlikely values of the random walk steps affect the validity of the construction 
of super- and sub-solutions described in Section [7} 

Remark 3. The error term in the statement of Theorem^ is 0(log(v* — v)), but 
we suspect that its correct order of magnitude might be 0(1), as suggested by the 
(mathematically non-rigorous) results obtained by B. Derrida and D. Simon in \YI\ 
[30] . and by analogous results obtained in the context of branching Brownian motion 
by J. Berestycki, N. Berestycki and J. Schweinsberg [6] and E. Aidekon, S. Harris, 
and R. Pemantle pp. Indeed, pushing the calculations just a little further with our 
method, it is not difficult to show that the error term is bounded above by a 0(1), 
i.e. that 

logPo(4») < + o(i). 

However, an extra log(v* — v) creeps into our lower bound, as a consequence of the 
sub-solution in Section^ being built with a(e) := — log(2) + e 2 instead of a(e) := 2, 
so that the best lower bound we can prove is 0(log(u* — v)). 
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6. Connection with the Brunet-Derrida theory of stochastic fronts 

Through a series of papers (see e.g. [T31 CEH EE ED EH 02]), Brunet and Derrida, 
partly with Mueller and Munier, developed a theory for the propagation of stochastic 
fronts described by perturbed F-KPP like equations. The methods and arguments 
used by Brunet and Derrida are not mathematically rigorous, but several of the 
corresponding results have now been given regular mathematical proofs (see e.g. 
[aSttgmESEHHHSlI!]). The purpose of this section is to explain how the 
present paper relates to these developments. 

We first explain how the survival probability of the branching random walk is 
related to travelling wave solutions of F-KPP like equations. We then explain how 
such travelling waves arise in the context of the Brunet-Derrida theory of stochastic 
fronts, and describe the approach used by Brunet and Derrida to deal with these 
travelling waves. We finally discuss how Theorem [21 its proof, and some of its con- 
sequencs, fit into the series of rigorous mathematical results establishing predictions 
from the Brunet-Derrida theory. 

6.1. Survival probability and travelling waves. The F-KPP equation, named 
after Fisher [20] and Kolmogorov, Petrovsky and Piscounov [26], is one of the clas- 
sical PDE models of front propagation, whose salient feature is to lead to travelling 
wave solutions. In its simplest form, the equation reads 

du d 2 u 

where u = u(x,t), x £ R, t > 0, and a travelling wave solution means a solution u 
of the form 

u(x,t) = g{x — vt), 
where «£lis the wave speed, and g(x), x£l describes the wave shape. 

In our context, it turns out that the survival probability qoo (as a function of the 
starting point x of the branching random walk) can be viewed as the shape of a 
travelling wave solution to a discrete analog of the F-KPP equation, with a special 
boundary condition. Indeed, define, for all n > and x£K, 

u n {x) := qoo(-x + nv). 

From Proposition [H one sees that the following equation holds for u n : for all x < 

(n + l)v, 

(29) u n+1 {x) = iJj{E(u n (x -())). 

Using the definition of ip, and rearranging the formulas a little, the above equation 
rewrite^ 

(30) u n+1 (x) - u n {x) = E{u n {x - C,)) - u n {x) + E{u n {x - ()) - E{u n {x - C)) 2 • 



dt lix 2 



Remember that the above equation holds only for x < (n + l)v, while, due to the 
fact that t/oo = on ] — oo,0[, we have that u n+ i(x) = for x > (n + \)v. It is 



^Written below the equation is the term-by-term analogy with the F-KPP equation. 
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now apparent that (n, x) h-> u n (x) describes a discrete travelling wave with speed 
v, obeying a discrete analog of the F-KPP equation at the left of the front, and 
identically equal to zero at the right of the front, where at time n the location of 
the front is nv. 

6.2. Brunet-Derrida theory of stochastic fronts. We begin by describing two 
distinct models considered by Brunet and Derrida. One is the stochastic F-KPP 
equation 



where W is a standard space-time white-noise, and N is large, with the initial 
condition u(x,t = 0) = 1]- O o,o]( a; )- 

The other is a particle system on the real line where a population of N particles 
evolves according to repeated steps of branching and selection. Branching steps 
are identical to those of the branching random walk considered in this paper, each 
particle in a given generation being replaced in the next generation by two new 
particles whose locations are shifted using two independent random walk steps. 
Selection steps consist in keeping only the ./V rightmost particles among the 2N 
obtained by branching from the population in the previous generation. 

It is possible to see these two models as describing the propagation of a front, 
and Brunet and Derrida found (see [13] [HI [H] ) that, for both models, the limiting 
velocity vn of the front has the following behavior as N goes to infinity: 

(32) Voo - U7V ~C(logiV)- 2 , 

where v^, is the limiting value of vn as N goes to infinity. Mathematical proofs 
of these results were then obtained in [27] for the stochastic F-KPP equation case, 
and in [4j for the branching-selection particle system case. Let us mention that the 
description of stochastic fronts obtained by Brunet, Derrida, Mueller and Munier 
goes far beyond ()32p . which is, in some sense, a first-order result (we refer to [111 
[101 [12] for more details). 

A very rough sketch of the argument used by Brunet and Derrida to deduce 
(l32j) is as follows. For solutions u of the stochastic F-KPP equation (}3"I) . at every 
time t > 0, x \— > u(x, t) continuously connects 1 at x = — oo to at some random 
x = X(t) defining the position of the front, right of which u(-,t) is identically zero. 
Looking at the equation (l3Tj) . one can see that stochastic effects due to the noise 
term counterbalance the u(l — u) creation term when u is of order 1/N. To find 
the asymptotic speed of propagation of the front, one should thus look for travelling 
waves obeying the F-KPP equation at the left of the front, taking values of order 
1/N near the front, and which are identically equal to zero at the right of the 
front. To study these travelling waves, one replaces the F-KPP equation by a linear 
approximation, for which explicit solutions can be found - these solutions should be 
approximately valid for the original equation, thanks to the fact that the values of 
u are small near the front. One can then check that the speed of these travelling 
waves must satisfy (I32p . The same line of argument is used to derive (I32p for the 
branching-selection particle system. In this case, the iV — > +oo limit of the model 
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can be viewed as a time-discrete version of the F-KPP equation similar to (|30p . and 
the actual system with N particles is described by a perturbation of this equation. 
As in the case of the stochastic F-KPP equation, one finds that the relevant scale 
for perturbations is 1/N, but this time this is due to the fact that a population of 
iV particles has a resolution of 1/N for representing a probability mass. 

Now, the key observation is that the dicrete travelling wave (n,x) t— >■ u n (x) de- 
duced from goo in the previous section corresponds just to the object studied by 
Brunet and Derrida in their argument for (132j) . namely travelling waves obeying an 
F-KPP like equation at the left of the front, that are identically equal to zero at the 
right of the front. One slight difference is that Brunet and Derrida prescribe the 
order of magnitude of the values of the travelling wave near the front to be 1/iV, 
and ask for the corresponding speed, while we prescribe the speed to be v = v * — e, 
and ask for the corresponding order of magnitude of the values near the front. 

Looking for a speed v(N) such that, for large N, 

Po(Ax>(A0)) ~ VN, 
we see from Theorem [T] that one must have 
(33) v* - v{N) ~ ^i*A"(t*)(log A0~ 2 , 

which corresponds precisely to the behavior (|32|) of vm obtained by Brunet and 
Derrida in the branching-selection particle system case, and is also the key to the 
rigorous proof of (|32p given in [3] . 

The connection between the survival probability of the branching random walk 
and the travelling wave solutions of perturbed discrete F-KPP equations of the type 
investigated by Brunet and Derrida, was in fact used by Derrida and Simon in 
|17[ [30] to derive Theorems [T] and in a mathematically non-rigorous way, from the 
Brunet and Derrida approach outlined above. 

6.3. Discussion. Our strategy for proving Theorem [2] is based on the original ar- 
gument of Brunet and Derrida sketched in the previous section, and turns out to 
be quite different from the probabilistic approach used by Gantert, Hu and Shi in 
|21j . Remember that the idea is to replace the discrete F-KPP equation by a linear 
approximation of it for which explicit solutions can be computed. In our context, 
this corresponds to replacing the non-linear convolution equation ([5]) by the linear 
equation (J7J). An additional idea we use is that the monotonicity properties of the 
non-linear equation allow for a rigorous comparison between suitably adjusted so- 
lutions of the linear equation, and solutions of the original non-linear one, via the 
construction of super- and sub-solutions. We were inspired by the work of Mueller, 
Mytnik and Quastel |27[ [28] , where this comparison idea is used as an intermediate 
step in the rigorous proof of (I32j) in the stochastic F-KPP equation case. In their 
context, the equation corresponding to (J5]) is the non-linear second-order ordinary 
differential equation — vv! = u" + u(l — u), and one can rely on specific techniques, 
such as phase-plane analysis, to implement the comparison idea. However, these 
tools are not available in our discrete-time setting, and we had to find a quite dif- 
ferent way of achieving the comparison argument. Note that, as a direct by-product 
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of the argument in [27\ [28] , one can prove an analog to Theorem [2] for the sur- 
vival probability of the branching Brownian motion killed below a linearly moving 
boundary. 

The asymptotic behavior of the survival probability qoo plays a key role in our 
proof of (|32p in the branching-selection particle system case, which was given in 
[4], and relied on the proof of Theorem Q] given by Gantert, Hu and Shi in |21| . 
In this regard, an interesting feature of the proof of Theorem [2] presented here is 
that, combined to the comparison argument described in [4], it provides a proof 
of (]32p in the branching-selection particle system case which is along the lines of 
the original argument by Brunet and Derrida. What is (a little) more, the slight 
improvement from Theorem Q] to Theorem [2] concerning the order of magnitude of 
the error term allows us to refine (|32p and thus achieve a result comparable with 
the one obtained in \27\ 128] . Indeed, under the assumptions of Theorem [21 one has 
that, in the branching-selection particle system case, 

v * _ VN _ ^*A"(O(logA0- 2 = O (^2$) • 

Note that the j^r§w term in the above equation corresponds to the actual order of 
magnitude expected from [IT], where more precise predictions are given. 
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Appendix A. The t* assumption 

Remember that the core results in this paper are derived under the following two 
assumptions on the random variable £ describing the steps of the branching random 
walk: 

(1) the probability distribution of ( has bounded support; 

(2) there exists t* > such that A(t*) - i*A'(t*) = -log(2), where A is the 
log-Laplace transform of £. 

The goal of the present section is to discuss the role and meaning of the second 
assumption, which we call the t* assumption. For the sake of simplicity, we keep 
assuming throughout this section that ( is bounded. We first try to provide some 
intuition about the t* assumption by giving some alternative characterizations. We 
then discuss the asymptotic behavior of the survival probability when the t* assump- 
tion is not met. 

A.l. Alternative characterizations of the t* assumption. The results con- 
tained in this section seem to be more or less folklore in the branching process liter- 
ature. However, we failed at finding a reference providing both complete statements 
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and self-contained proofs, so we decided to include the following short elementary 
account of these results. 

We start by giving a general characterization of the critical speed v*. Let £+ be 
the essential supremum of £, and let £_ be its essential infimum. Then let 9 : R — > R 
be defined by 

9{t) = log (2£(exp(iC))) = log(2) + A(t). 

From the bounded support assumption for £, 9 is well-defined and finite for all t € R, 
and is C°° and convex as a function of the parameter t. Then define 

(34) v* = mi!^^,t>0 

From results by J. M. Hammersley [22J, J. F. C. Kingman |25j, and J. D. Big- 
gins [7], v* appears as the maximum limiting speed in the branching random walk. 
Indeed, letting M n denote the position of the rightmost particle of the branching 
random walk at time n, one has that 

Theorem 3. With probability one, — > v* . 

We refer to [8] for a survey of this type of results. Note that, as an immediate 
consequence, £_ < v* < C+- 

The following result shows that the existence of a t* such that ([2]) is satisfied 
(the t* assumption), is equivalent to the infimum defining v* in (|34p being in fact a 
minimum. (Note that (|36p below is exactly the same as ([2j).) 

Lemma 4. For all t* > 0, 

(35) e Jn =v . 

is equivalent to 

(36) ^ = 9>(t*), i.e. Kit*) - t*A'(t*) = - log(2). 

Proof. Let t* > 0. If ([35]) holds, then the map defined by t \-t 6(t)/t reaches its 
minimum at t*. But its derivative at t* is : 

9'(t*)t* -9{t*) 
t* 2 

Therefore (I36p holds. Conversely, assume that ([36]) holds. By convexity of 9 and by 
(I36|) we get, for all t > 0, 

e(t) > e'{t*)(t - f) + 9(t*) = 9(t*)t/t*. 

Therefore (ESI) holds. □ 



The following lemma provide some alternative probabilistic interpretations for the 
t* assumption. 

Lemma 5. The following conditions are equivalent. 

(1) There exists t* > such that Ii36\) holds. 

(2) P(( > v*) > 0, i.e. v* < C+. 
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(3) P(C = C+) < 1/2- 

Proof. By Lemma [TJ Condition 1 implies Condition 2. Let us check that Condition 
2 implies condition 1. We have: 

lim — — = +oo and lim — — = C+- 

Assume v* < f + . Then, v* can be defined as the infimum of t \— > 9(t)/t on a compact 
interval of ]0, +oo[. Therefore, there exists t* > such that (|35p holds. By Lemma 
0]we then get that Condition 1 holds. 

Let us check that Conditions 2 and 3 are equivalent. This is a special case of 
Proposition A. 2 in [23]. Below we give a sligthly modified proof. Let us write 

d(t) = t(+ + logCA + B(t)) 

where 

A = 2P(C = C+) and B(t) = 2E(exp(t(( - C+))k<C + )- 
Note that A + B(t) converges toward A. If P(( = ( + ) < 1/2, we then get: 

lim log(A + £(£)) < 

(the limit might be — oo). Therefore there exists t > such that 8(t)/t < £+ and 
we then have v* < £+: Condition 2 holds. If, on the contrary, P(£ = £ + ) > 1/2, 
then for all t > we have A + > 1 and then 6(t)/t > C+. Therefore v* > ( + : 
Condition 2 does not hold. □ 

A. 2. Survival probability when the t* assumption is not fulfilled. Our aim in 

this subsection is to show that, when the t* assumption is not fulfilled, the behavior 
of the survival probability differs from the one established in the core of the paper. 

Let us consider the process of children displaced by exactly C+ from their parent 
in the branching random walk. This is a Galton- Watson process. The number of 
children of a given parent follows a binomial distribution with parameters 2 and 
P(£ = The behavior of the surving probability ¥q(A oq (v)) as v tends to v* 

from below turns out to depend on the behavior of this Galton- Watson process. 

(1) Subcritical case, i.e. P(£ = C+) < 1/2. By Lemma[5]this is the case studied 
in the paper. We have: 

logPo(^ooW)~-C C (^-u)- l/2 

where Cq is an explicit positive constant that depends on the distribution of 

C- 

(2) Critical case, i.e. P(( = C+) = 1/2. By Lemma [5] we have v* = (+. As in 
the subcritical case we have: 

P (4x») ^P (^oo («*)), 

since Po(^4 o( , w*)) = is also the survival probability of the critical Galton- 
Watson process. However, the rate of convergence is not the same as in the 
subcritical case. Indeed, there exists a positive constant depending only 
on the distribution of ( such that 

n(Aoo(v)) >d c (v*-v). 
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When P(£ = 0) = P(£ = 1) = 1/2 this is the lower bound given by Propo- 
sition 2.4 of [29]. The general case follows by coupling and rescaling@. 

The precise behavior of the survival probability in fact depends on the 
behavior of P[v < ( < v*). When there exists v < v * such that P(v < £ < 
v*) = there exists a constant such that 

n(Aoo(v))<D c (v*-v). 

When P(( = 0) = P(( = 1) = 1/2 this is the upper bound given by 
Proposition 2.4 of [29|. We get the previous result by coupling and rescaling 
as above. 

When, on the contrary, P(v < £ < v*) > for all v < v * the survival prob- 
ability can be larger. Indeed, set p v = P(C > v) > 1/2. The Galton- Watson 
process obtained by keeping all steps greater or equal to v is supercritical. 
Denote by q v > its survival probability. We have 

(1 -p v +p v (l - q v )) 2 = (1 - q v ). 

Therefore 

_ 2p v - 1 

Qv — q 
Pv 

We then have: 

2r> — 1 

Po(Ax») > qv = 8P(v<C<v*) 

pi 

(let us recall that P(( = v*) = 1/2 and that P(( > v*) = 0). Thus, the 
survival probability Pq(A oq (v)) can tend arbitrarily slowly to 0. 
(3) Supercritical case, i.e. P(C = C+) > 1/2- By Lemma[5]we have v* = £+, as 
in the critical case. Let us denote by q > the survival probability of the 
Galton- Watson tree. Then, we have: 

Po(A>») -»• Po(AooK)) = ^ > 0. 
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First replace £ by a random variable £ such that P(£ = ^+) = P(£ = f_) = 1/2, thus lowering 
the survival probability. Then replace £ by (£ — C-)/(C+ — C-)- 



